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ABSTRACT 

Two methods were used to calculate the heat flux 
to full-coverage film cooled airfoils and, subsequently, 
the airfoil wall temperatures. The calculated wall 
temperatures were compared to measured temperatures 
obtained in the Hot Section Facility operating at real 
engine conditions. Gas temperatures and pressures up 
to 1900 K and 18 atm with a Reynolds number up to 
1.9 million were investigated. Heat flux was calcu- 
lated by the convective heat transfer coefficient adia- 
batic wall method and by the superposition method which 
incorporates the film injection effects in the heat 
transfer coefficient. The results of the comparison 
indicate the first method can predict the experimental 
data reasonably well . However, superposition overpre- 
dicted the heat flux to the airfoil without a signifi- 
cant modification of the turbulent Prandtl number. 

The results of this research suggests that additional 
research is required to model the physics of full- 
coverage film cooling where there is significant tem- 
perature/density differences between the gas and 
coolant. 

NOMENCLATURE 

hg heat transfer coefficient defined by Eq. (1) 

hg ( e ) heat transfer coefficient defined by Eq. (5) 

K. constant 

L total surface length from leading edge 

stagnation point 



p pressure 

q heat flux 

Re Reynolds number 


s equivalent slot height 

T temperature 

V/V cr critical velocity ratio. 



x surface distance from stagnation point 

Xl distance from applicable film cooling row 

y specific heat ratio 

n film cooling effectiveness 

e temperature difference ratio, Eq. (6) 

u viscosity 

p density 

Subscripts 

aw adiabatic wall 

c coolant out 

9 gas 

ge effecti ve gas 

meas measured 

w wal 1 outer surface 

1,2,3, rows of film cooling holes 
INTRODUCTION 

Improved performance of turbojet and turbofan 
engines is typically accompanied by increased cycle 


pressure ratio and combustor exit gas temperature. Gas 
pressure levels of 25 to 30 atm and gas temperatures of 
1600 K exist in some current operational engines while 
pressure levels up to 40 atm with gas temperatures of 
1800 K are anticipated in advanced commercial engines. 
These continuing increases in the turbine entry gas 
pressure and temperature of the modern gas turbine 
engine and the high development cost puts a premium on 
an accurate initial aerothermal design of the turbine 
hot section hardware. 

The design goals for commercial jet engines 
include high cycle efficiency, increased durability of 
the hot section components (lower maintenance costs), 
and lower operating costs. These goals are contradic- 
tory in that high cycle efficiency reguires minimizing 
the cooling air reguirements while increased durability 
reguires metal temperatures and temperature gradients 
to be minimized. An optimum design can only be real- 
ized through an improved understanding of the flow 
field and the heat transfer process in the turbine gas 
path. 

Sophisticated computer design codes are being 
developed which have the potential of providing the 
designer with significantly better initial estimates 
of the flow field and heat load on the hot section com- 
ponents. These codes are being evaluated and verified 
through low temperature and pressure research in cas- 
cades and tunnels. However, by design, these facili- 
ties do not model all of the processes that exist in a 
real engine environment, and therefore, the ability of 
the design codes to predict the interaction of the 
various parameters cannot be fully evaluated. 

The Hot Section Facility at the NASA Lewis 
Research Center provides a "real-engine" environment 
with established boundary conditions and convenient 
access for advanced instrumentation to study the aero- 
thermal performance of turbine hot section components. 
The thermal performance and, ultimately, the life of 
these components in a realistic application is depen- 
dent on the designer's ability to predict the local 
heat load distribution. 

The heat flux to a stator airfoil, which has full 
coverage film cooling, presents a particularly chal- 
lenging situation for the designer because of the com- 
plex flow field and the interaction between the hot and 
cold air streams. Even though film cooling has been 
studied for several years, there is a limited amount 
of realistic engine data available to verify the ana- 
lytical models. 

Two models have been developed to predict the heat 
flux to film cooled surfaces and have been discussed 
and compared in the literature (1) . A widely used 
method defines the heat flux based on standard convec- 
tive heat transfer coefficient and an adiabatic wall 
temperature. The influence of film cooling is then 
described in an effectiveness parameter. This approach 
is discussed in Ref. 2. For full-coverage film cool- 
ing, another approach is to redefine the heat transfer 
coefficient by incorporating the effects of coolant 
injection in the coefficient and then computing heat 
flux using the recovery temperature of the gas stream 
{3 ) . In almost all verif ication/eval utation experi- 
ments used to define the correlation coefficients for 
these methods, low temperature, constant property, 
heated surface models are considered. In real engine 
applications, these deviations can have a significant 
effect on the designer's ability to predict the temper- 
ature and life of the structure. 

The purpose of this paper is to examine these two 
approaches to calculating heat flux at high temperature 
and compare predicted wall temperatures of a turbine 
airfoil with measurements taken in the Hot Section 
Facility. These experimental data were taken at 


Reynolds numbers of 1.25 and 1.90 million, gas tempera- 
tures and pressures of 1350 to 1900 K and 9 to 18 atm, 
and coolant blowing ratios from 0.5 to 1.5. 

The results are presented as a comparison of ana- 
lytical and experimental airfoil temperature distribu- 
tions. Each method of calculating the heat flux to a 
full coverage film cooled surface is compared to the 
experimental data for each Reynolds number conditions. 


General description 

A physical layout of the Hot Section Facility 
(HSF) is shown in the perspective view (Fig. 1(a)). 

The HSF facility located at NASA Lewis Research Center 
is a unigue facility having fully-automated control of 
the research rig through an integrated system of mini- 
computers and programmable controllers. The major 
components of this facility and how they interface 
together to provide a real engine environment are shown 
in the flow diagram (Fig. 1(b)). This facility is dis- 
cussed in more detail in Refs. 4. and 5. 

The main air supply system provides air at 10 atm 
to a nonvitiated preheater. The preheater modulates 
the air temperature between ambient and 560 K. Through 
a set of routing valves, two modes of operation can be 
selected. Utililzing the compressor bypass system, air 
can be provided to the test rig at 10 atm and up to 
560 K. The second mode, compressor-mode, can provide 
air to the research rig at pressure up to 20 atm and a 
temperature up to 730 K when utilizing the heat of 
compression. 

The research test rigs (Fig. 1(a)) consist of two 
independent test stands: a modified turbine test rig 

(full annular cascade rig) and a combustor test rig. 

The combustor test rig was used to develop and document 
the exit temperature profile and efficiency of the heat 
source for the full annular cascade. Utilizing three 
pairs of instrumentation rakes (temperature, pressure, 
and exhaust products) located at the exit of the com- 
bustor, the circumferential and radial profiles were 
documented providing a known input profile for the 
cascade vane row. 

Cascade configuration 

A cross-section of the Hot Section Cascade Rig is 
shown in Fig. 2. The major components consist of a 
heat source (combustor), the full annular vane row 
(containing full-coverage film cooled vanes), an 
exhaust duct line, a guench system (to lower the tem- 
perature of the exhaust gas), and the exhaust system. 

The vane row consists of 36 full-coverage film- 
cooled (FCFC) stator vanes. The 36 vanes are separated 
into two groups: 10 tests vanes and 26 slave vanes. 

The test vane and slave vane cooling air is supplied 
from two separate manifolds with the flow rates to each 
manifold independently computer controlled. 

FCFC vane 

The stator vane configuration used for these tests 
was a full-coverage film cooled design with an impinge- * 

ment insert to provide augmented coolant-side heat 
transfer. The vane row hub and tip diameters were 
0.432 and 0.508 m, respectively. Both the vane height 
and chord were 3.81 cm. More detailed geometric data w 

are given in Table I and Ref. 6. 

A typical slave vane and test vane are shown in 
Fig. 3(a). The cooling air supply tube on the tip of 
the test vane allows cooling air to be supplied from a 
manifold separate from the slave vane supply. The test 
vane shown in Fig. 3(a) is in its finished form while 
the slave vane is shown in a partially finished form. 
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A detailed airflow calibration of each test vane 
was performed before installation in the stator case. 
This information was used to select the ten test vanes 
with similar flow characteristics and to provide 
"pressure loss" coefficients for the computer code. 

The stator case is shown partially assembled in 
Fig. 3(b). The ten test vanes and some instrumentation 
leads are clearly shown. The cavity directly over the 
vane row feeds cooling air to the slave vanes while a 
separate manifold (downstream of the stator row) feeds 
cooling air to the test vanes. The top dead center 
positon (zero degree) is also noted in the figure. All 
circumferential locations arc measured in a counter 
clockwise (CCW) direction from this position lookinq 
downstream. 

Instrumentation 

Research instrumentation in the Cascade Rig con- 
sists primarily of the conventional steady-state 
temperature and pressure measurements. The vane air- 
foils carry thermocouples and pressure sensing tubes 
to sense the gas side metal temperatures, gas-stream 
static pressures, and cooling air side temperatures and 
pressures. Cooling air flow is controlled and measured 
at the venturi in each of the supply lines. Cooling- 
air temperatures and pressures are also measured in the 
internal manifolds of each cooling air system. 

Gas path conditions are monitored at the stator 
row inlet and exit. Each location has three fixed 
probe ports for mounting radially-actuated water-cooled 
probes. Radial gas-path surveys of both temperature 
and pressure are recorded from vane hub to tip. 

A cross-sectional schematic of the vane airfoil is 
shown in Fig. 4 with a composite summary of instrument 
locations on the airfoil. The locations shown in 
Fig. 4 represents either metal temperature or static 
pressure measurements. Because each airfoil could 
accommodate only a limited number of instrumentation 
grooves, the temperature or pressure distributions 
reported are composed of measurements from several 
airfoils in the test vane sector of the annulus. 

EXPERIMENTAL PROCEDURE 

There were two basic modes of facility operation: 
10- and 20-atm. Within each mode the research riq was 
operated with the combustor operating (burning) or 
without combustor operation (isothermal). The research 
objectives were to investigate the aerothermal perform- 
ance of this cascade over a range of Reynolds numbers, 
and at a constant Reynolds number, over a range of com- 
bustor exit temperatures and pressures. The various 
operating modes of this facility is shown in Fig. 5 and 
Table II. The data reported herein are taken from 
cases 10 and 12 . 

The gas conditions were established by setting the 
combustor inlet total pressure, the vane exit outer 
radius static pressure, and the combustor fuel/air 
ratio through predetermined input values stored in the 
operations computer. The coolant flow rate and temper- 
ature were varied systematically at fixed gas condi- 
tions either through predetermined input values to the 
computer or by manually inputting values to the 
computer. 

ANALYTICAL PROCEDURE 

The two methods that were used to determine heat 
flux to the turbine airfoil will be presented first. 
Then the method of calculating the airfoil wall temper- 
atures will be discussed. 


Gas-side heat flux 

Method I. TTTi s method defines the heat flux using 
a conventional convective heat transfer coefficient and 
an adiabatic wall temperature. 

q = h (T - T ) (1) 

H g ' aw w 

where the adiabatic wall temperature is defined by a 
film effectiveness equation. 



where Tg e is the effective gas temperature based on 
recovery a of the velocity. 

Although it is recognized that injection into the 
boundary layer will affect the heat transfer coeffi- 
cient, particularly near the injection site, the value 
of this coefficient is typically equated to the 
"unblown" value. The effects of injection, geometry, 
etc. are incorporated in the film effectiveness term. 

The STAN5 (_7) boundary layer code was used herein to 
compute the "unblown" heat transfer coefficient on the 
airfoil. The boundary layer was assumed to be "tripped" 
to a turbulent state at the first row of holes. 

The film effectiveness was based on a "typical" 
correlation format found in the literature. One form 
of correlation that was developed for convex-curved 
surfaces is (2) 



Even though this form was originally developed for 
slots it is also used for rows of holes. The constant 
K was found to be 1.50 in Ref. 8 which investigated 
two rows of film coolant injection onto a convex-curved 
surface with test conditions that simulated engine 
operation. 

For a full-coverage film cooled surface, the film 
effectiveness related to each succedinq row of holes 
in the downsteam direction is influenced by the film 
from upstream rows. The film effectiveness for the 
furthest downstream row can be calculated by the 
following summation (9) • 

ifilm = n i + n 2 ^ ~ + ^3 (1 _ n^)(l - n 2 ) (A) 

where nfiim is the composite effectiveness down 
stream of the last row of holes. 

Method 2 - The effect of film cooling is incorpo- 
rated in the heat transfer coefficient for this method 
and the heat flux is defined by the following equation 
derived in Ref. 3. 

i - y> tv - v (5) 

where a dimensionless temperature parameter is 
defined as 



Through linear superposition arguments the value 
of the heat transfer coefficient at any value of e 
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Eq. 3. This approach resulted in wall temperatures 
calculated for the suction surface which have a slope 
similar to the experimental temperatures but with a 
magnitude 100 to 200 K higher. 

Kasagi et. al. (_13) presents results that show the 
film effectiveness on the pressure surface to be lower 
than that of the suction surface for a given "blowing" 
ratio. Consequently, the calculations for the pressure 
surface were made with the constant of Eq. (3) set at 
1.20. The calculated wall temperatures show good 
agreement with the experiment both in slope and 
magnitude. 

The results of method 2 calculation are also shown 
in Fig. 9. Similar trends as method 1 are noted. The 
pressure surface calculations compare with the experi- 
ment reasonably well on the forward portion of the air- 
foil but deviate up to 100 K higher on the aft portion. 
The calculations on the suction surface follow the same 
pattern as method 1. In order to obtain this compari- 
son it was necessary to increase the turbulent Prandtl 
number by a factor of 3.5 above the freestream value 
of 0.86 as the boundary layer calculation approached 
the wall. A factor of 2.0 is recommended in Ref. 7, 
however, a very poor comparison between analysis and 
experiment was obtained using 2.0. Variations in DELMR 
and ALAM of 2:1 from the recommended values were unable 
to improve the comparison. The two empirical param- 
eters, DELMAR and ALAM, which are incorporated into the 
STANCOOL model to account for film injection were not 
sufficient to model the real engine situation. 

There is no clear-cut reason why the turbulent 
Prandlt number should increase by this much as it 
approaches the wall. However, the experimental data 
in the literature used to determine its functional 
relation are generally measured with a "hot" wall and 
"cold" air stream which tends to destabilize the boun- 
dary layer. The data of the experiment reported herein 
were taken on a "cold" wall with a "hot" gas stream. 

This high rate of cooling of the gas stream is usually 
associated with a stabilizing effect on the boundary 
layer. This is consistant with the lower heat transfer 
coefficients calculated with the higher turbulent 
Prandtl number near the wall. This phenomena does not 
apply to the method 1 calculation since the correlation 
(Eq. (3)) is based on experiment with a "cold wall" and 
a "gas hot" stream. 

There has been some arbitrariness in selecting 
empirical coefficients for the calculation procedures. 
The primary guide for the selection has been the exper- 
imental data. However, there is a basis for the selec- 
tion. The point to be made is that neither calculation 
method has modeled the physics well enough to permit an 
accurate prediction of airfoil temperatures without 
some prior knowledge. Method 1 does model the suction 
surface data reasonably well especally in trend which 
is not surprising since the correlation is based on 
convex curvature data. Method 2, which is more rigor- 
ous in accounting for the physics of the problem, does 
as well as method 1 on the suction surface once the 
turbulent Prandtl number is modified. However, 
method 2 does not predict with any degree of accuracy 
the phenomena on the airfoil aft pressure surface. 

There are many additional parameters in the 
STANCOOL model which could affect the calculated heat 
flux to the airfoil surface. Even the assumption of 
linearity of the energy equation could be investigated. 
However, the purpose of this paper is to point out the 
state of predicting temperatures on a full-coverage 
film-cooled airfoil in a realistic, hostile environment. 
In addition, topics for additional research with 
properly scaled temperature/density ratios have become 
apparent. 


SUMMARY OF RESULTS 

Two methods were used to calculate the heat flux 
to full-coverage film cooled airfoils and, subse- 
quently, the airfoil wall temperatures. The calculated 
temperatures were compared to measured temperatures 
obtained in the High Pressure Facility. The following 
results were obtained. 

The convective heat transfer coefficient, adia- 
batic wall temperature approach predicted the slope of 
the experimental airfoil temperatures reasonably well 
on the suction surface. However, the magnitude of the 
prediction was up to 200 K higher than the experimental 
data. The pressure surface temperature data was pre- 
dicted reasonably well by decreasing the empirical 
coefficient of the film effectiveness correlation by 
20 percent. 

The modified superposition approach did as well 
as the first method on the suction surface for both 
slope and magnitude. However, there was a relatively 
poor comparison with the experimental data on the aft 
portion of the pressure surface. 

The superposition method required a modification 
of the turbulent Prandtl number near the wall in order 
to obtain the relatively good comparison with the 
experimental data and with method 1. The results of 
this research suggests that additional research is 
required to model the physics of full-coverage film- 
cooling with significant temperature/density differ- 
ences between the gas and coolant as found in a gas 
turbine environment. 
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TABLE I. - STATOR VANE GEOMETRY 


Mean diameter, cm 46.99 

Vane height, cm 3.81 

Axial chord, cm 3.81 

Axial solidity 0.929 

Aspect ratio 1.000 

Number of vanes 36 

Leading edge radius, cm 0.508 

Trailing edge radius, cm 0.089 


TABLE II. - HSF CASCADE RESEARCH CONDITIONS 


Case 

Combuster exit 


Exit 

Mode of 
operation , 
atm 



■ 

■ 

Reynolds 

number 




Isothermal 

mode 






1 

500 

5.1 

0 . 

23 

320 

0.775 

1 . 25x10® 

10 

2 

500 

6.5 







1.90xl0 6 

10 

3 

500 

8.5 







2.50xl0 6 

10 

4 

700 

4.9 







1.25xl0 6 

20 

5 

700 

7.4 







1.90xl0 6 

20 

6 

700 

9.5 







2.50xl0 6 

20 





Burning 

mode 







1300 


0 . 

23 

' ■ 

0.775 

0.5x10® 

||H| 


1535 




■ 

.775 

0.5xl0 6 



945 




• 

.71-. 83 

1.90x10® 

KJ»W m 

mm 

1165 





.72-. 80 

1.25x10® 


■ 

Bl 

1500 





.775 

1.25x10® 



HI 

1500 





.775 

1.90x10® 
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(b) HSF flow schematic. 
Figure 1. - Hot section test facility. 
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(a) The slave vane is shown in a partially machined form while the 
test vane is shown in finished form. 

Figure 3. - Full coverage film cooled vanes for the Hot Section Facil- 
ity cascade rig. 
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(b) Partial assembly of FCFC stator case showing test vane sector. 

Figure 3. - Concluded. 
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Figure 4. - Full coverage film cooled vane cross-section showing typical 
composite of instrumentation. 



COMBUSTOR 

EXIT 

PRESSURE, 

atm 


DESIGN OPERATING 



COMBUSTOR EXIT TEMPERATURE, K 


REYNOLDS NO. IS BASED ON VANE EXIT CONDITIONS AND THE VANE CHORD. O/A 
IS THE BOUNDARY LAYER COOLING RATE C °"’ 


Figure 5. - Cascade rig simulation of real engine operating conditions. 


CRITICAL VELOCITY RATIO, V/V 



Figure 6. - Critical velocity ratio design conditions for the 
stator airfoil tested. 
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(a) Suction surface. (b) Pressure surface. 

Figure 8. - Film-cooled heat transfer coefficients for method 2 calcu' 
lated by the Stancool boundary layer code. 
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(a) Reynolds number * 1 250 000. (b) Reynolds number - 1 900 000. 

Figure 9. - Comparison of airfoil temperatures calculated by method 1 and method 2 with experimental data. 
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